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SECTION  I 


INTRODUCTION 


The  United  States  Air  Force  requires  the  capability  to  operate  aircraft 
on  rough  runways  created  by  the  rapid  repair  of  bomb  damage  (Reference  l). 
The  current  rapid  runway  repair  technique  uses  AM-2  matting  and  procedures 
described  in  Air  Force  Regulation  93-2  (Reference  2).  Initial  data  for  F-4E 
performance  on  AM-2  mats  have  been  acquired  in  the  HAVE  BOUNCE  test  program 
(References  3  and  4).  These  data  can  also  be  used  to  validate  computer 
simulations  of  aircraft  response  to  runway  surface  roughness.  These  com¬ 
puter  simulations  may  then  be  used  to  develop  surface  roughness  criteria 
(Reference  5 ) . 

The  United  States  Air  Force  Academy  (USAFA)  has  provided  engineering 
consultation  on  surface  roughness  to  the  Air  Force  Engineering  and  Services 
Center  (AFESC)  since  1977.  During  this  reporting  period,  evaluations  have 
been  provided  for  using  statistical  analysis  techniques  and  using  a 
stationary  test  facility  for  surface  roughness  evaluations.  The  purpose 
of  this  report  is  to  develop  a  three-degree-of-freedom  model  for  analysis 
of  aircraft  taxi  response  over  repair  mats  and  spalls.  A  one-degree-of- 
freedom  model  was  used  to  predict  the  symmetric  motion  response  for  the  HAVE 
BOUNCE  Phase  II  program  (Appendix  C,  Reference  4). 

This  report  uses  modal  analysis  to  predict  aircraft  response  to  surface 
roughness.  The  rigid  body  aircraft  model  includes  the  bounce,  pitch,  and 
roll  degrees  of  freedom.  The  surface  roughness  for  the  nose  and  two  main 
gear  tracks  is  modeled  as  a  sequence  of  straight  line  segments,  such  as 
ramps  and  mats,  placed  on  an  originally  flat  runway  surface.  The  equations 
of  motion  are  integrated  using  modal,  superposition  and  the  piecewise  exact 
solution  algorithm  (Reference  6).  The  model  will  accept  experimentally 
determined  natural  frequencies,  damping  ratios,  and  mode  shapes.  The  modal 
superposition  solution  is  suitable  for  programming  on  minicomputers.  The 
results  presented  include  a  comparison  with  the  TAXI  computer  code  (Refer¬ 
ence  7)  and  predictions  of  asymmetric  motion  due  to  spalls  in  the  runway 
surface . 
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SECTION  II 


EQUATIONS  OF  MOTION 


1*  GEOMETRY.  The  aircraft  model  geometry  and  coordinates  are  shown 
in  Figure  1.  The  spring  and  damping  constants  represent  the  strut  and  tire 
combined  in  series  for  vibrations  about  the  equilibrium  position.  The 
attachment  point  displacements  are  kinematically  related  to  the  center  of 
gravity  displacements  through  a  geometry  matrix: 


Figure  1.  Aircraft  Geometry  and  Coordinates 
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2.  EQUATIONS  OF  MOTION.  The  equations  of  motion  are  derived  assuming 
nose  profile  and  attachment  point  displacements  are  ordered  so  that 


>  u,T  >  0 


“HP  >  “N 


and  their  slopes  are  ordered  so  that 


>  u„  >  0. 


"  UP 


Then  the  total  upward  nose  strut  force  acting  on  the  aircraft  is 


FSN  ^N  ^NP  “  “n^  +  CN  ^P  "  ^N^ ' 


Similarly,  the  upward  main  gear  strut  forces  are 


FSR  *R  ^P  "  ^  +  cr  Wp  " 

FSL  =  *1,  ^ULP  '  +  CL  ^ULP  ’  ^ 


The  equations  of  motion  from  summing  forces  vertically  and  summing 
moments  about  the  pitch  and  roll  axes  are 


fsn  +  fsr  +  fsl  ~  mog  uog 


dnf sn  '  dm  ^fsr  +  ip  9 


Define 


D  F  -  D  F  =  I 
L  SL  R  SR  R 


MCG  0  0 


[m]  =  0  Ip  0 


0  0  Ir 


These  equations  can  be  written  in  matrix  form  as 


Cm]{^}  +  =  C^lCupJ  +  [c^Cup}  (U) 


where 


PN 


=  <  UPL  ) 


u_ 

PR 


[K_]  =  [GEOM]] 


*N  ° 

0  *1 
0  0 


0 

0 

*R. 


[C  ]  =  [GEM*]1 


CN  0 


0  cT 


0  0 


0 

0 


“RJ 


Then  using  Equation  (D  the  left-hand  side  of  Equation (4)  can  be  re¬ 
written  in  center-of -gravity  displacements  as 


M(^)  -  CCpltvy.)  +  tKp1( Uj,}  -  [Kjltu,,}  + 


(5) 


where 


[Kp]  =  [GEOM]1 


[Cp]  =  Cgeom]1 


*N 

0 


N 

0 

0 


*L 

0 

0 

cT 


0 

0 

^R 

0 

0 


[GEOM] 


0  cT 


[GEOM]  . 


vi  '<»•  ,-v. 


3.  MODAL  EQUATIONS.  The  modal  equations  of  motion  are  derived  using 
the  natural  frequencies  and  mode  shapes  for  the  undamped  system  which  are 
computed  from  standard  vibrations  techniques  as  in  Reference  8  or  computer 
techniques  as  in  Reference  9.  The  system  characteristic  matrix  equations 


(Dtj,]  -  4  m)  f  MODEi}  =  {0} 


where  go*  is  the  natural  frequency  and  { MODE^}  is  the  corresponding  mode 
shape.  Let  [MODAL]  be  the  modal  matrix  constructed  using  the  mode  shapes 
as  columns.  The  mode  shapes  are  orthonormalized  so  that 

[modal]t[m][modal]  *  [I]  (7) 

where  [i]  is  the  identity  matrix.  Then 

[modal]t[kd][modal]  =  3  (8) 

is  a  diagonal  matrix  with  the  square  of  the  natural  frequencies  on  the 
diagonals . 

The  same  sequence  of  matrix  products  does  not  generally  diagonalize 
the  damping  matrix 

[modal]t[cd][modal]  =  [Cp]  (9) 

It  is  very  desirable  that  a  diagonal  matrix  [Cp]  be  obtained  so  that  the 
equations  of  motion  will  decouple  with  the  introduction  of  principal  co¬ 
ordinates  {p}  defined  by 

[uG]  =  [MODAL] [p]  (10) 

For  lightly  damped  systems  such  as  seen  in  the  HAVE  BOUNCE  Phase  I 
test  results  (Reference  3)>  it  is  reasonable  to  neglect  the  coupling  of 
modes  through  the  off-diagonal  terms  of  [Cp]  and  retain  only  the  diagonal 
terms  (Reference  6).  These  diagonal  terms  are  then  rewritten  as 

Vi =  2  9  “i 

where  is  the  ratio  of  damping  in  mode  i  to  the  critical  damping  for  that 
mode. 

The  decoupled  modal  equations  of  motion  are  obtained  through  multi¬ 
plying  Equation(5)by  [MODAL]^  and  substituting  Equation  (10),resulting  in  the 
sequence  of  equations  of  the  form 


■  «* ;  .«?4* 


(12) 


p±  +  2  ?i  WjPi  +  W^Pi  =  {M0DEi3T([K1]{up}  +  CC^fUp}) 

Thus,  the  response  of  each  mode  can  be  solved  as  a  single  degree  of  freedom 
system  and  the  total  response  found  by  superposition  of  the  principal  co¬ 
ordinate  responses  using  Equation (10) . 
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SECTION  III 


SOLUTIONS 


1.  DIRECT  INTEGRATION.  There  are  several  methods  for  directly 
integrating  the  linear  matrix  equations  of  motion  [Equation  53  such  as  the 
Newmark  Method  and  Wilson  0  Method  (References  6  and  9).  Some  of  these 
methods  are  conditionally  stable, and  each  requires  a  small  time  step  to 
minimize  period  elongation  and  amplitude  decay  of  the  results. 

2.  MODAL  SUPERPOSITION.  This  method  integrates  the  decoupled 
equations  of  motion  [Equation (12 )]  separately  and  then  adds  each  modal 
response  for  the  total  solution.  If  a  direct  integration  method  such  as 
the  Wilson  0  Method  with  the  same  time  step  for  all  modes  is  used,  the  re¬ 
sults  will  be  the  same  as  directly  integrating  the  system  matrix  equations. 
But  modal  superposition  provides  the  opportunity  to  only  integrate  the 
modal  response  for  selected  dominant  modes  and  to  use  the  accurate  piece- 
wise  exact  integration  scheme  which  does  not  require  a  small  time  step  for 
stability  and  accuracy  as  do  the  direct  Integration  schemes. 

3.  PIECEWISE  EXACT  SOLUTION  METHOD.  Many  types  of  loading, including 
the  runway  mat  profile, can  be  closely  approximated  as  a  sequence  of 
straight-line  segments  with  unequal  time  intervals.  For  example,  Figure  2 
shows  a  runway  mat  profile  with  upheaval  and  settlement  as  constructed  for 
the  HAVE  BOUNCE  Phase  II  tests  (Reference  4).  Each  straight-line  segment, 
such  as  the  second  up-ramp  illustrated  in  Figure  2,  has  its  elevation 
written  as  a  linear  equation  of  the  form  e, (t)  =  a^  +  bjt  where  the  sub¬ 
script  i  denotes  the  segment  number.  The  time  axis  values  t  are  the  product 
of  multiplying  the  runway  distance  x  and  the  constant  ground  speed  of  the 
aircraft  V. 


Figure  2.  Piecewise  Linear  Runway  Profile 

For  any  segment  of  time  ts  <  t  <  t^,  the  one -dimensional  vibration 
equation  to  be  integrated  is 

x  +  2  C  COX  +  w2x  =  CO2  (a  +  bt)  =  A  +  Bt  (13) 

where  a  +  bt  is  the  equation  for  the  mat  profile  in  this  time  interval.  The 
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exact  displacement  solution  to  Equation  (13) is  (Reference  6) 


x(t)  =  Aq  +  +  A^e  ^  cos  co^t  +  A^e  ^  ^  sin  w^t 


where 


A  2  Sb 


A  =  —  - 

0  2  .3 

0)  W 


Ai=S 

0) 


a2  =  *(V  -  A0 

A3  =k,  (i(t0}  +  C  wA2  -  Ax) 
d 

Wd  =  uji  -  C2 


(14) 


The  solution  is  exact  so  using  a  small  time  step  for  stability  and 
accuracy  is  not  required. 

4.  ONE-DIMENSIONAL  MODAL  ANALYSIS.  The  HAVE  BOUNCE  Phase  I  test 
results  indicate  the  aircraft  responds  in  a  dominant  bounce  mode  which  is 
lightly  damped  (Reference  3).  Based  on  these  results,  it  appears  that  a 
single  degree  of  freedom  solution  based  on  Equation(l4)would  give  good 
estimates  for  main  gear  compression  loads.  The  one-dimensional  equation 
of  motion  is 

u  +  2  £  u)  u  +u)2u  =  u)2uM(t)  +  2  C  tou^t)  (15) 

where  u«(t)  is  the  main  gear  track  profile.  The  frequency  U3 and  damping 
factor  C  may  be  obtained  experimentally  from  a  power  spectral  density 
analysis.  The  second  term  on  the  right-hand  side  of  Equation (15) may  be 
neglected  due  to  the  small  amount  of  damping  in  the  system.  The  main  gear 
compression  load  is  the  sum  of  the  dynamic  forces  and  the  weight  (W). 

Fm  =  M^j  u)2(taM(t)  -  u(t))  +  Mg,  2  C  a)  ((^(t)  "  u(t))  +  w  (16) 

For  lightly  damped  systems,  the  second  term  may  be  neglected  and  the  re¬ 
sulting  Fm  will  be  a  conservative  estimate  of  the  gear  load. 

This  model  has  been  programmed  at  the  Air  Force  Flight  Test  Center 
(AFFTC)  and  shows  good  correlation  with  the  HAVE  BOUNCE  Phase  II  test  re¬ 
sults  (Reference  4). 


5.  THREE-DIMENSIONAL  MODAL  ANALYSIS.  The  success  of  the  one¬ 
dimensional  model  indicates  that  a  three-dimensional  model  could  provide 
good  results  for  the  bounce,  pitch,  and  roll  degrees  of  freedom.  This 
model  would  estimate  nose  gear  forces  as  well  as  compute  asymmetric  roll 
response. 

The  modal  analysis  is  based  on  integrating  Equation  (13) using  Equation 
(l1*) for  each  mode  and  each  gear  track  profile.  For  three  modes  and  three 
tracks,  this  solution  will  require  integrating  Equation (14) up  to  nine 
separate  times  and  adding  the  results. 

Since  Equation  (14)  is  being  integrated  in  principal  coordinates,  the 
gear  loads  are  not  entered  directly.  Instead,  the  gear  loads  are  parti¬ 
tioned  using  the  mode  shape  vector.  This  partitioning  accounts  for  how 
much  the  gear  load  is  exciting  a  particular  mode  shape. 


For  example,  consider  the  modal  solution  due  to  left  main  gear  profile 
only.  The  gear  load  is  partitioned  using  the  orthonormali zed  mode  shapes 
that  satisfy  Equation (7).  Hie  first  modal  equation  is 


px(t)  +  2  Cj^Pj/t)  +  o^P^t)  =  PLl(t) 
where  the  load  participation  is 

PL-Jt)  =  {mode1}t[k1]  ^.(t)} 

0 


(17) 


(18) 


For  each  segment  of  time  where  pL^(t)  =  A  +  Bt,  the  solution  for 
pn(t)  is  found  using  Equation  (l4)  with  co  =  u>j.  The  load  participation  for 
the  other  modes  is  found  using  the  appropriate  mode  shape  in  Equation (l8) 
and  the  response  found  using  the  appropriate  frequency  and  damping  factor 
in  Equation  (l4i  The  center-of -gravity  displacements  are  found  using  the 
modal  matrix 


UCG 

pi 

1  6 

)  =  [MODAL]  ( 

P2  > 

w 

! 

and  the  attachment  point  displacements 


(19) 


are  found  using  the  geometry  matrix 


H 

r  -\ 

UCG 

k 

}  =  [GEOM]  { 

9  > 

w 

0 

The  gear  compression  loads  are  computed  using 
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?L  =  (UL  "  V  +  h 


FR  =  -kE  \  +  FMST 


where  Fjjg<j>  and  F^gy  are  nose  and  main  gear  static  loads.  The  response  for 
the  nose  and  right  main  tracks  follow  in  a  similar  manner. 
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SECTION  IV 


RESULTS 


Calculational  results  for  the  single-degree -of -freedom  model  (SDOF) 
and  the  three-degree-of-freedom  model  (3D0F)  will  be  compared  with  results 
from  the  TAXI  computer  code  (Reference  7).  Then  spall  results  using  the 
3D0F  model  will  be  presented. 

1.  TAXI  CODE.  This  code  was  provided  by  the  Air  Force  Wright 
Aeronautical  Laboratory  (AFWAL).  Additional  FORTRAN  coding  was  added  to 
interface  with  the  USAFA  computer  plotting  system  and  to  evaluate  slope  of 
the  strut  air  curve  formulas  to  obtain  spring  constants  for  the  linear 
vibrations  model. 

2.  COMPARISON  MODEL.  The  heavy-weight  aircraft  data  of  Table  1  was 
obtained  from  AFFTC.  The  initial  sprinq  constants  for  the 

vibrations  models  were  obtained  from  the  TAXI  code  by  computing  the  slope 
of  the  strut  air  curve  at  the  static  equilibrium  position.  This  slope 
gives  an  overall  stiffness  for  the  strut  which  is  then  added  to  the  tire 
stiffness  in  series  (Reference  8)  to  obtain  an  overall  gear  stiffness. 


TABLE  1.  AIRCRAFT  MODEL  DATA 


mcg 

=  145.5  lb-sec2/in 

Ip 

=  1.92(106)  lb-sec2* in 

XR 

=  6.56(105)  ib-sec2-in 

°L  "  °R 

=  109-25  in 

dn 

=  23^. 9  in 

dm 

=  44.35  in 

TABLE  2.  EQUIVALENT  SPRING  STIFFNESS 


Main  Gear 

Strut  Stiffness 

27809  ib/in 

Tire  Stiffness 

13530  lb/in 

Series  Stiffness 

9102  lb/in 

Nose  Gear 

Strut  Stiffness 

1662  lb/in 

Tire  Stiffness 

15600  lb/in 

Series  Stiffness 

1502  lb/in 
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3.  T11REE-DEGREE-OF-FREEDOM  ( 3DOF)  MODEL .  The  linear  vibrations 
model  for  the  data  of  Tables  1  and  2  has  the  mode  shapes  and  frequen¬ 
cies  shown  in  Table  3. 


TABLE  3.  MODE  SHAPES  AND  FREQUENCIES 


Mode  1 

Mode  2 

Mode  3 

co  7.271  rad/s 

12.016  rad/s 

18.198  rad/s 

xg  0.025931 

-0.078743 

0 

8  0.000685 

0.000226 

0 

0  0 

0 

0.001235 

The  mode  shapes  indicate  relative  motion  within  a  mode  and  are  orthonormalized 
according  to  Equation  (7)  .  The  first  mode  has  its  node  37.9  inches  aft  of  the 
CG  so  its  primary  motion  is  pitching  at  the  nose  gear.  The  second  mode  has 
its  node  348.4  inches  forward  of  the  CG  so  its  primary  motion  is  bouncing. 

The  third  mode  is  the  roll  mode  which  is  decoupled  frctn  the  bounce  and  pitch 
modes  due  to  symmetry  of  the  aircraft. 

4.  SINGLE -DEGREE -OF -FREEDOM  (SDOF)  MODEL.  This  model  is  based  on 
the  bounce  mode  so  it  uses  U3 =  12.016  rad/ s  and  M  =  145.5  sec2/ in.  Both 
SDOF  and  3DOF  models  were  programmed  using  CAL  (References  6  and  10) ,  a  com¬ 
puter  analysis  language  developed  for  educational  use  at. the  University  of 
California  at  Berkeley.  The  basic  programming  of  Equation  (14)  is  simple 
and  can  be  done  in  almost  any  computer  language. 

5.  RESPONSE  TO  HALF  MAT  PROFILE.  The  half  mat  profile  used  for  the 
following  calculations  is  shown  in  Figure  3.  The  TAXI  code  response  for 

a  constant  velocity  of  20  knots  is  shown  in  Figure  4.  The  frequency  of 
vibration  after  coming  off  the  mat  is  1.80  Hertz, and  the  ratio  of  damping 
to  critical  damping  is  about  1%. 

The  SDOF  response  for  a  natural  bounce  mode  frequency  of  1.91  Hertz 
(Mode  2)  and  damping  factor  of  1  percent  is  shown  in  Figure  5.  The  straight 
line  through  the  plot  is  the  static  equilibrium  load.  The  maximum  compres¬ 
sion  load  is  sliqhtly  larger  than  the  TAXI  response,  but  the  minimum  values 
are  much  smaller  including  unloading  at  the  fourth  minimum. 

The  3DOF  response  for  a  20-knot  constant  run  with  1  percent  damping 
in  all  three  modes  is  shown  in  Figure  6.  The  straight  lines  through  the 
plots  are  static  equilibrium  values.  The  maximum  main  gear  compression 
loads  are  slightly  smaller  than  both  the  TAXI  and  SDOF  results,  and  the 
minimum  values  lie  in  between  the  other  two  runs.  The  nose  gear  response 
is  smaller  to  the  TAXI  response  but  does  not  contain  the  highest  frequency 
motion. 


6.  RESPONSE  TO  FULL  MAT  PROFILE.  The  full  mat  profile  used  for  the 
following  calculations  is  shewn  in  Figure  7.  The  TAXI  code  response  for  a 
20-knot  run  is  shown  in  Figure  8.  The  main  gear  compression  load  coming 
off  the  mat  (at  three  seconds  on  the  figure)  shows  a  reduction  of  the  am¬ 
plitude  of  vibration.  A  comparison  of  the  TAXI  code  responses  of  Figures  4 
and  8  indicates  that  the  trailing  edge  of  the  mat  can  cause  reinforcement 
or  cancellation  of  the  motion.  To  obtain  accurate  simulations  of  experi¬ 
mental  results  showing  such  reinforcements  and  cancellations,  the  computer 
model  must  have  the  same  period  of  vibration  in  its  dominant  mode  as  the 
aircraft . 

The  SDOF  response  for  a  frequency  of  I.91  Hertz  and  damping  factor  of 
1  percent  in  shown  in  Figure  9.  The  initial  response  coming  on  the  mat  is 
similar  to  the  TAXI  response,  but  the  linearized  frequency  of  1.91  Hertz 
causes  the  phasing  of  the  trailing  edge  of  the  mat  so  that  the  motion  is 
ampiif iea. 

The  3DOF  response  for  a  20-knot  run  with  1  percent  damping  in  all  three 
modes  is  shown  in  Figure  10.  The  main  gear  response  is  very  similar  to  the 
SDOF  response  showing  an  amplification  at  the  end  of  the  mat.  The  nose  gear 
response  is  very  similar  to  that  for  the  TAXI  code  response. 

The  difference  in  phasing  at  the  end  of  the  mat  is  due  to  a  different 
effective  frequency  of  vibration  about  the  equilibrium  position.  The  SDOF 
and  3D0F  results  used  the  slope  of  the  air  curve  at  equilibrium  as  a  com¬ 
ponent  value  to  compute  the  series  stiffness  of  the  gear  assembly.  The 
TAXI  code  uses  a  nonlinear  model  for  the  compression  air  curve  behavior  of 
the  strut  as  well  as  nonlinear  damping. 


Figure  3-  Half  Mat  Profile 

To  evaluate  the  effects  of  the  nonlinear  TAXI  model,  the  effective 
spring  constant,  the  change  in  pneumatic  force  divided  by  the  change  in 
strut  compression  was  computed  for  each  output  value.  The  main  gear  spring 
constant  ranges  from  a  low  of  U733  lb/in  to  a  high  of  11,362  lb/in.  The 
nose  gear  spring  constant  ranges  from  a  low  of  965  lb/in  to  a  high  of  ?6b? 
lb/in.  This  variation  in  spring  constants  should  be  averaged  in  such  a  way 
that  the  equivalent  spring  constant  yields  the  same  period  of  free  vibration 
while  the  aircraft  is  on  the  mat.  Since  the  period  is  inversely  proportional 
to  the  square  root  of  the  spring  constant,  it  is  reasonable  to  compute  the 
equivalent  spring  constant  according  to 


r 

I 


The  revised  spring  constants  are  8081  lb/in. for  the  main  gear  and  1508  lb/in 
for  the  nose  gear.  The  revised  mode  shapes  and  frequencies  are  shown  in 
Table  4. 


TABLE  4.  REVISED  MODE  SHAPES  AND  FREQUENCIES 


Mode  1 

Mode  2 

Mode  3 

U) 

7.282 

11.328 

17.148 

xg 

0.025052 

-0.079027 

0 

0 

0.000688 

0.000218 

0 

0 

0 

0 

0.001235 

The  bounce  mode  frequency  is  now  1.80  Hertz.  The  SDOF  results  using 
this  fundamental  frequency  are  shown  in  Figure  11.  The  3D0F  results  using 
the  revised  mode  shapes  and  frequencies  of  Table  4  are  shown  in  Figure  12. 
Figures  11  and  12  now  show  much  better  agreement  with  the  TAXI  code  results 
of  Figure  8.  in  particular,  the  cancellation  of  motion  at  the  trailing 
edge  of  the  mat  is  consistent  for  each  calculation. 


7.  SPALL  MODEL  RESULTS.  The  spall  profiles  used  for  the  following 
calculations  are  shown  in  Figure  13.  These  profiles  correspond  to  the 
spalls  cut  into  the  runway  for  testing  at  AFITC.  The  spalls  were  cut  across 
the  runway  so  that  tests  could  be  made  with  the  left  main  gear  only  entering 
the  spall  and  other  tests  could  be  made  with  both  the  nose  gear  and  left 
main  gear  entering  the  same  spall  profile. 

Figure  14  shows  the  response  for  the  left  main  gear  entering  the  two- 
foot  spall  with  a  speed  of  10  knots.  Figure  15  shows  the  response  for  the 
nose  gear  and  left  main  gear  entering  the  two-foot  spall  cut  across  the 
runway  with  a  speed  of  10  knots.  Figure  16  shows  the  symmetrical  response 
for  all  three  gears  entering  the  two-foot  spall.  The  main  gear  compression 
loads  are  of  about  the  same  magnitude  as  shown  for  the  mat  encounters.  The 
maximum  left  main  gear  compression  load  is  larger  for  the  asymmetric  spall 
input  (Figure  15)  than  for  the  symmetric  spall  input  (Figure  16). 

Figure  17  shows  the  response  for  the  left  main  gear  entering  the 
double  spall  with  a  speed  of  10  knots.  Figure  18  shows  the  response  for 
the  nose  gear  and  left  main  gear  entering  the  double  spall  cut  across  the 
runway.  Figure  19  shows  the  response  for  all  three  gears  entering  the 
double  spall.  The  sharp  downward  plunge  of  the  main  gear  compression 
follows  from  the  assumption  that  the  gear  follows  the  runway  profile.  The 
effect  of  the  tire  footprint  following  such  a  spall  profile  is  not  in¬ 
cluded  in  these  results  but  may  be  added  by  modifying  the  runway  profile 
input.  The  maximum  main  gear  compression  load  is  larger  for  the  asymmetric 
input  (Figure  18)  than  for  the  symmetric  input  (Figure  19). 
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Figure  13.  Spall  Models 
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SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  one-  and  three-degree-of-freedom  linear  vibration  models  developed 
in  this  report  can  reproduce  the  response  of  the  essentially  nonlinear  air¬ 
craft  system  to  surface  roughness.  To  correctly  reproduce  the  response  to 
runway  mats,  it  is  important  that  the  linear  model  have  spring  constants 
that  produce  the  same  period  of  free  vibration  so  that  the  phasing  of  the 
trailing  edge  of  the  mat  will  cause  the  same  reinforcement  and  cancellation 

phenomena . 

The  three-degree-of-freedom  model  can  accept  arbitrary  runway  profile 
inputs  for  each  gear.  It  can  be  used  to  predict  results  of  asymmetric 
motion  due  to  spalls  in  the  runway.  The  predictions  of  this  model  need  to 
be  validated  using  the  results  of  the  recent  spall  tests  at  the  AFFTC. 
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